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Abstract. We discuss a cosmological model in which the string gauge field coupled univer- 
sally to matter gives rise to an extra centripetal force and will have effects on cosmological 
and astronomical observations. Several tests are performed using data including galaxy ro- 
tation curves of twenty- two spiral galaxies of varied luminosities and sizes, and perihelion 
precessions of planets in the solar system. 

Remarkable fit of galaxy rotation curves is achieved using the one-parameter string 
model as contrasted to the three-parameter model of dark matter model with the Navarro- 
Frenk- White profile. The rotation curves of the same group of galaxies are independently 
fit using dark matter model with the generalized Navarro- Frenk- White profile and using the 
string model. The average of the NFW fit is 9% better than that of the string model at a 
price of two more free parameters. From the string model we give a dynamical explanation 
of Tully-Fisher relation. We are able to derive a relation between field strength, galaxy size 
and luminosity, which can be verified with data of the 22 galaxies. 

To test the hypothesis of the universal existence of string gauge field, the string model is 
also used to explain the anomalous perihelion precession of planets in the solar system. The 
required field strength is deduced from the data. The field distribution resembles a dipole 
field originated from the Sun. A refined model assumes that the string gauge field is partly 
contributed by the Sun and is partly contributed by the stellar matter in the Milky Way. 
The field strength needed to explain the excess precession is of similar magnitudes as the 
field strength needed to sustain the rotational speed of the sun inside the Milky Way. 
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Introduction 

Recent cosmological and astronomical observations are becoming increasingly interesting 
laboratories for precision tests of new physical theories aiming at extending the standard 
paradigms. On the one hand high energy completion of gravity theory should leave signatures 
sufficiently different from the low energy effective theories the like of general relativity and 
its modifications. These will be detectable with the advances in detector technologies. On 
the other hand, string theory and other quantum gravity candidate theories when applied to 
cosmology or astronomy, should shed new light on old problems. 

In this work we will confine our interest to two problems. The first one is the "missing 
mass problem" in galaxies-the mass discrepancy between the one measured by rotational 
speeds of stars inside a spiral galaxy and the one predicted from its stellar matter distribution. 
Another one is the anomalous precession of the planets inside the solar system recently 
reported, the explanation of which cannot be found within the standard framework. The 
model we propose to solve these two problems at the same time is a very special string 
model. The model, first discovered by Nappi and Witten [1], falls into a general class of 
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exactly solvable string models but it has the added merit that all effects due to the finite 
size of the strings are taken into account. In other words it is a bona fide string model. 
And it lives in four dimension spacetime which resembles as closely as possible to Minkowski 
spacetime but with the presence of the string gauge field. Due to the existence of the string 
gauge field, the geodesies in this four dimensional space-time are concentric circles instead 
of the usual straight lines in Minkowski spacetime. This is analogous to the Laudau orbits 
of an electron in the presence of a magnetic field. Thus Cheung et al [2] proposed the model 
to explain the galaxy rotation curves in spiral galaxies in lieu of Cold Dark Matter (CDM). 
In this paper we are going to extend their work with a direct comparison of goodness of fit 
of the string model with the CDM model in galaxy rotation curves fitting. Furthermore we 
are going to apply the idea to explain the anomalous precession of planets' perihelions in the 
solar system. 

1 Galaxy rotation curve 

While the Cold Dark Matter cosmology has been accepted by many as the correct theory 
for structure formation at large scale and the solution to the missing mass problem at the 
galactic scale there still lacks hard proof for the existence of Dark Matter particles. Until 
the coming Linear Hadron Collider and future experiments-see exciting development in this 
direction [3, 4]-tells us definitely what constitute Dark Matter a more natural and universal 
explanation in lieu of dark matter cannot be excluded. Here we entertain the possibility that 
a higher-rank gauge field universally coupled to strings can give rise to a Lorentz force in 
four dimensions providing an extra centripetal acceleration for matter towards the center of 
a galaxy in addition to the gravitational attraction due to stellar matter. If not properly 
accounted for, it would appear as if there were extra invisible matter in a galaxy. A salient 
feature of this Lorentz-like force is that it fits galaxies with an extended region of linearly 
rising rotational velocity significantly better than the dark matter model. This feature also 
endows the model with testability: in the region where the gravitational attraction of the 
visible matter completely gives way to the linear rising Lorentz force, typically in the region 
r ~ 20Rd, should we still observe linear rising rotation velocity, say for the satellites of the 
host galaxy, it would be a strong support for the model. Otherwise if all rotation curves are 
found to fall off beyond 20i?rf for all galaxies then the model is proven wrong. Furthermore 
this is the first of such attempts to verify directly the validity of string theory as a description 
of low energy physics. Given this last reason alone we regard it a worthwhile endeavour. 

1.1 The string model 

Consider a four-dimensional string model proposed by Nappi and Witten [1] in which the 
string theory is exactly integrable. Furthermore tree-level correlation functions-describing an 
arbitrary number of interacting particles-have been computed, which capture all finite-sized 
effect of the strings to all orders in string scale [5]. This is valid for all energy scales as long 
as the string coupling constant is weak. This is the region of interest when we extrapolate 
to our low energy world. The three-form background gauge field, -ff(3), coupled uniquely to 
the worldsheet of the strings, is constant. 

Because we are no longer approximating strings as point particles, this coupling between 
the two- form gauge potential B^2) the two dimensional worldsheet of the string produces 
a net force on the string when it is viewed as a point particle [5]. The center of mass of the 
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closed string executes Landau orbits: 



a = ao + re^^* (1.1) 

where a is the complex coordinate of the plane in which the time-like part of the three-form 
field has non-zero components. This phenomenon is analogous to the behaviour of an electron 
in a constant magnetic field. 

For the purpose of data fitting we only need to retain a component of the tensor gauge 
field which is perpendicular to the galactic plane of the spiral galaxies. We denote the 
strength of this gauge field by H, which together with the "charge-mass" ratio, forms the 
only free parameter of this model, denoted by = ^H. This is to be contrasted with the 
three free parameters one needs in the dark matter model using the celebrated Navarro- 
Frenk- White profile [6] for an iso-thermal and isotropic dark matter distribution, (the 
characteristic radius of the dark matter distribution), ctq (the central density), and a (the 
steepness parameter). A general profile with a free parameter a is used because of the need 
to fit rotation curves from dwarf galaxies as well as galaxies with a varied surface brightness, 
from high surface brightness to low surface brightness. 

The simplistic property of the string model is, in fact, a favoured approach from the 
string theory point of view, because the coupling of the field to matter has universal strength, 
i.e. all matter is charged, rather than neutral. Therefore if all matter is indeed made up 
of fundamental strings and hence couples universally to the tensor gauge field, H, each star 
in a galaxy will then execute the circular motion in concentric landau orbits on the galactic 
plane. Effectively there is an additional Lorentz force term in the equation of motion for a 
test star: 

m — = qHv+mF^ (1-2) 
r 

where the field, H, is generated by the rotating stellar matter and the halo of gases alike. 
1.2 The Dark Matter Model 

According to the Cold Dark Matter (CDM) paradigm there is approximately ten times more 
dark matter than visible matter. The fiuctuations of the primordial density perturbations 
of the universe get amplified by gravitational instabilities. Hierarchical clustering models 
furthermore predict that dark matter density traces the density of the universe at the time 
of collapse and thus all dark matter halos have similar density. Baryons then fall into the 
gravitational potential created by the clusters of dark matter particles, forming the visible 
part of the galaxies. In a galaxy the dark matter exists in a spherical halo engulfing all of 
the visible matter and extending much further beyond the stellar disc. To describe the dark 
matter component we use the generalized Navarro-Frenk- White profile: 

~ (JL^ad _|_ Jl)3-a ' 

where a = 1 corresponds to the NFW profile, and is the characteristic radius of the dark 
matter halo. In the Dark Matter fitting routine we allow the to vary from 3Rd to 30Rd- 
We further require that the dark matter density be strictly smaller than the visible mass 
density, ao < pQ. (The data can in fact be fit equally well with the roles of dark matter and 
visible matter inverted.) Here we treat do, Vg and a as free parameters. Together with po 



-3- 



and Rd from the visible component, the Dark Matter model utilizes five free parameters. All 
in all the rotation velocity of a test star is given by 



- = F, + Fdm (1.4) 
r 

in the Dark Matter model. 
1.3 Stellar matter 

To describe the visible matter we use the parametric distribution with exponential fall off in 
density due to van der Kruit and Searle in both models: 

p{r,z) = poe sech?{^) (1.5) 

with pq being the central matter density, Rd the characteristic radius of the stellar disc and 
Zd the characteristic thickness. Following a common practice we choose Zd to be \Rd] the 
dependence of the final results on this choice is very weak [7] . 

Gravitational attraction due to the visible matter is henceforth given by 

F,{r) = GNPoRdF{f) (1.6) 

where after rescaling f = 

Fir) = I f --'^--hyz') 



df 



all space 



becomes a universal function for all galaxies. 

Summary of equations in both models: We are now ready to fit the galaxy rotation 
curves data. In the string model the three free parameters, il, Rd, and po are defined by the 
following equation: 

— = GNpoRdF + 2^v . (1.7) 
r 

The fundamental charge-to-mass ratio and the strength of the gauge field is encoded alto- 
gether in one free parameter = f^- 

The five free parameters in CDM model are defined by 



GNPoRdF + FoM . {U 



The two free parameter po ^ind Rd are common in both models, describing stellar contribution 
to the rotational velocities of stars about the center of galaxy. FjjMii") is given by the 
following expression: 

FDM{r) =4ttJ dr'- 



(^)°(1 + ^) 

All in all the CDM profile carries another three free parameters, namely, the central density of 
dark matter halo, ctq, the characteristic length scale of the halo, Rg, as well as the "steepness" 
parameter of the halo, a. 



-4- 



1.4 Fitting Procedures: 

A few remarks concerning our fitting procedures are in order. The Dark Matter model with 
its five free parameters and the string model with its three free parameters are independently 
fit to the data to obtain the best fit values for each set of the parameters for each of the 
twenty-two spiral galaxies. Under no circumstances, the best fit values from one model are 
fed into the other model as prior values. Except restricting the to vary from 3Rd to 30Rd 
in CDM model and letting to be gi?^ in the stellar distribution as conventionally done 
(see for example [7]) to save computing time there are no other simplifications. We then use 
the best fit values for these two parameters (and three others in the dark matter model) to 
compute the total mass, as well as the mass-to-light ratios, for these galaxies. These will 
serve as sanity check for the best fit values of the parameters in both models. 

Independent of any galaxy rotation modelling, p and Rd can be fit with photometric 
data and hence are not really free parameters. Since we are only interested in comparing the 
dark matter model and our string model in fitting the galaxy rotation curves we are treating 
them as free parameters in each model. We instead choose to use our best fit values for these 
parameters, from Dark Matter model and String model in turn, to compute the total mass in 
each galaxy and cross check with independent astronomical observations. We also compute 
the percentages of baryonic matter in the galaxies for the CDM model as usually done by 
astronomers. This serves as another check of our methodology. 

We are ignoring the gas contributions from our fitting; because putting in more free 
parameters will no doubt improve the fit for both models. For the same reason we do not 
allow for any correction for star extinctions and supernova feedback as they would not affect 
any conclusion we draw concerning the relative quality of the fit between the CDM model and 
string model. Keeping this simplistic spirit we do not allow for any dark matter component 
at all in the string model and we also assume that the strength of the string field be constant 
throughout the span of each galaxy. Local back reaction of spacetime to the presence of the 
string field is also ignored. 

1.5 Analysis 

The rotation velocity of a given test star is solved from equations (1.7) and (1.8) for the string 
and the dark matter model, respectively. The data are then fit to the dark matter model 
and string model independently. The best fit values of the parameters from each model are 
obtained by minimizing the functionals. ^We obtained our rotation curve data of the 
twenty-two galaxies in the SINGS sample from the FaNTOmM website. Using the best fit 
values of these parameter of the dark matter model, we can compute the mass of the dark 
matter halo, the mass of stellar mass, and then the ratio of luminosity to the total mass. 
These values are tabulated in Table C of Appendix C. From the string model we can likewise 
determine the set of values for the strength of the string gauge field, the mass of stellar mass 
and then finally the ratio of luminosity to the stellar mass. These values are tabulated in 
Table 9 in Appendix C. 

^Notice that in the observations the distance measurement from the center of the galaxy is assumed to be 
exactly. The uncertainty is attributed, instead, to the velocity measurements. During fitting, however, we 
discovered that uncertainty in determine the "center of the galaxy" affects the quality of the fit significantly. 
According to both models the rotation velocity at the "center" of the galaxy should be exactly zero. If we 
could shift some data by a linear translation to make the zero velocity point coincide with the r = point 
by hand, we would have obtained much lower values for both models. Therefore this linear shift is better 
attributed to the error in distance determination. 
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Figure 1. Rotation curve of NGC2403 fit with the dark niatter model (left) and with the string 
model (right). The x-squared value per degree of freedom using the dark matter model with a NFW 
profile is 4.515 while that using the string model is 4.304. The X-axis is radius in kpc and Y-axis is 
velocity in kms~^. 

In Figure 1 rotation curve of galaxy NGC2403 fitted using NFW profile (left) and 
string model (right) are plotted side by side for comparison. Squares with error bars are 
observational data. Theoretical predictions are indicated by the solid lines with stars in the 
NFW fit (left) and with triangles in the string fit (right). The string model clearly gives a 
better fit. 

The x-squared value, per degree of freedom, from the string fit is 4.304 whereas that 
from the NFW fit is 4.515. Overall string model gives a x-squared value of 1.656 averaged 
over the 22 galaxies while the dark matter model gives 1.594. The fitting results of 22 galaxies 
using dark matter model and string model are detailed in Appendix C. The best fit values 
of the free parameters are tabulated in Table C for Dark Matter model and in Table 9 for 
String model, respectively, in Appendix C. We can see that the NFW profile fits marginally 
better at a price of two more free parameters. 

After we obtain the best fit values for the free parameters we can compute the (total) 
masses for the galaxies. 

String Model: For this model there is only visible matter whose mass can be straightfor- 
wardly computed by integrating (1.5) with the best fit values of Rd and po for each galaxy. 

NFW profile: Matter in this model consists of the visible matter, same as that in the 
string model, and the dark matter which assumes the generalized NFW density profile (1.3). 
The NFW profile gives divergent mass if the radius is integrated to infinity. We therefore 
adopt the usual cutoff and compute the mass only up to the virial radius within which the 
average density is 200 times the critical density for closure. 

1.5.1 Visible Mass to Light vs B-magnitude 

Using the measured B-band absolute magnitudes we compute the visible mass to light ratios 
for the galaxies. In string model these ratios fall between 0.11 ~ 5.6 and centred around 1 
as shown in Fig. 3. The same ratios from the NFW model span five orders of magnitude 
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Figure 2. The total-mass-to-light ratios derived from the best fit values and the raeasured B-band 
luminosity of the 22 galaxies in dark matter model. 




Figure 3. The visible-mass-to-light ratios derived from the best fit values and the measured B-band 
luminosity of the 22 galaxies in the string model. 



(see Fig. 2) with a lot of them falling much below 1. For the NFW model we also compute 
the percentage of baryonic matter in the total mass. According to the CDM paradigm this 
number should be around 10%. However the actual results come with a wild scatter. The 
scatter in the mass-to-light ratios and the baryon fractions clearly indicate that NFW profile 
is not capturing the underlying physics correctly. 

1.5.2 TuUy-Fisher relation 

A Tully-Fisher relation can be derived from the string model which relates the rotation 
velocity in the "flat" region of the rotation curves to the product of the total luminous mass, 
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M^, and the parameter, Q . 



From the equation of motion (1.2) we solve for v, 



(1- 



:i.io) 



We then look for a balance of falling Newtonian attraction and rising Lorentz force, 
resulting in ^ ~ 0. Because we know that the turning point is at r ~ 2.2rd, setting ~ 
yields a relation between and 
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Inside the orbit r ~ 2.27?^ lies most of the visible mass. We can therefore use the point- 
mass approximation when computing and Upon substituting (1.11) into (1.10) our 
TuUy-Fisher relation follows. The string model therefore provides a dynamic origin of this 
well-tested rule of thumb. 
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Figure 4. The luminous mass and velocity relation of the 22 galaxies fit by the Tully-Fisher relation 
derived from the string model. 



We plot our best fit values of G M* Q against v in Fig. 1.5.2. The representative velocity, 
V, is selected to be the maximal observed velocity in the entire curve for each galaxy, to 
eliminate man-made bias. This no doubt introduces more scatter than necessary. Despite 
that the data obey the relation remarkably well. 

Note that this is a nontrivial relation for it relates two parameters from two additive 
force terms to an observed quantity, the rotational speed. Furthermore if one can determine 
the luminous mass of the galaxy, M^, and the field strength, H, we can determine a funda- 
mental ratio, the charge-to-mass ratio, universal for all matter according to string theory, by 
measuring the rotating speed of a test star. Furthermore our model provides a dynamical 
explanation to the Tully-Fisher relation. 
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1.5.3 A relation obtained from the string model fitting result 

By dimensional analysis, together with some common results from astronomy, we can find 
a simple relationship between the field strength Q, galaxy luminosity and size. This serves 
as a consistency check for the string model. According to the string model, considering its 
analogy with electromagnetism, it is reasonable to expect the average field strength to be 
proportional to M"i?~^ where M is the total luminous matter in the galaxy and R is the size 
scale of the galaxy. ^ To be more specific we will heavily use the electromagnetism analogy in 
the following discussion. Consider a group of electrons azimuthal symmetrically distributed 
and in rotation around the z axis. Let us look at the magnetic field at the center of this 
distribution, the determining physical quantities are: the magnetic constant fj,Q, mass density 
scale po, distribution size scale Rq and rotational angular velocity scale coq. ^ 

(Other determining factors include the shape and the spatial dependence of the mass 
distribution, the spacial distribution of the angular velocity. These factors do not change 
the result of dimensional analysis but change the proportional constant.) By dimensional 
analysis, we have 

B (X Ho ■ po ■ ujQ ■ RI (1.12) 

Using the total charge Q (x po ■ Rq, and defining vq = loq ■ Rq, we have 

Bocpo-Q-Ro^ -vo (1-13) 

Note that vq is the rotational velocity scale for the galaxy. Translating to the language of 
the string model, it is 

QckQ-Rq"^ -vo (1.14) 

The proportional constant here depends only on the mass and angular velocity distributions, 
or abstractly on the galaxy type. ^ Thus galaxies with similar mass distribution profile and 
rotation curve shapes should have similar constants of proportionality. Now recall M (x Q, 
where the proportional constant is universal and thus the same for all galaxies. Furthermore 
we also use the assumption L oc M. ^ Thus 

UocL-Rq^-vo (1.15) 

To relate vq to L we use the Tully-Fisher relation which says L oc A V"' where a is around 
2.5 lb 0.3 and AV is the velocity width of the galaxy[8]. The proportional constant in the 
Tully-Fisher relation is galaxy type independent. Since vq is the overall scale for AV, we 
have also L oc (fo)"- Using this in (1.15), we get 

n-Rlo^L^+^ (1.16) 

^By the Tully-Fisher relation [8] velocity is related to the total mass of the galaxy, therefore we do not 
need a separate term for the velocity dependence. 

^What we really want to check is the averaged field over the galaxy, but it is proportional to the field 
strength at the center. 

''More precisely, it is not just the galaxy morphology type. The velocity distribution also matters. It 
is possible galaxies of the same morphology type but very different velocity distributions will have different 
proportional constants. 

^This relation is independent of the galaxy type. For more discussion about the mass luminosity-relation 
among galaxies of different types, see [9]. 
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Figure 5. luminosity to ilR^ log-log plot 



which after taking logarithm 
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where k is the proportional constant in the relation (1.16). The log-log diagram is shown in 
Fig.5. 

There seems to be a trend of a linear relation in the "main" part of the diagram. The 
slope of the line is about |, quite close to the 1 -|- ^ 3 derived theoretically. Two points 
seems to lie outside of the "main" part, i.e one at the lower left corner for the dwarf galaxy 
mSldwb, another one for NGC4236 at the left of the upper right group. In terms of galaxy 
type these two galaxies are "exotic" among the 22, and thus perhaps their Inn deviated more 
from those in the "main" part. Actually in terms of morphology type, NGC4236 is SBdm, 
which is the most irregular one among the regular types, while mSldwb is the only dwarf 
galaxy among the 22 galaxies [10], all others are more or less regular galaxies. Presumably 
these two lie on another line for irregular galaxies which is parallel to the line passing the rest. 
It is possible that there are series of parallel lines for different types of galaxies. However, 
statistical error from the small size of this data set makes above arguments weak. Analysis 
of this kind for a larger number and more types of galaxies could make the situation clearer. 
But this exercise is beyond the scope of the current project. 

1.6 Discussion 

The original appeal of the NFW profile based on the ideas of hierarchical clustering was its 
universality. One simple NFW profile was expected to explain structure formation, rotation 
curves of galaxies, giant or dwarf, from high to low surface brightness. This promise has 
been undermined by the cusp and core debate in dwarf galaxies as well as in the low surface 
brightness galaxies (see for example [7]). The fact that light does not follow dark matter-well 
established by detailed observation and analysis in the Milky Way (see [11] for a summary), 
in addition to a clear deficit of satellite galaxies in MW have only served to thicken the 
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plot. Recently this debate has been taken to broader context by observational progress: The 
simplicity of the galaxies [12] and the early assemble of the most massive galaxies [13] are 
at odd with the hierarchical clustering paradigm. While we are not claiming that our string 
toy model can answer all these questions in one stroke we merely show that it pans out just 
as well as the Dark Matter model in fitting the galaxy rotation curves while using two fewer 
parameters. Moreover by tuning the ratio of the strength of the string field to stellar mass 
density galaxies with a wide range of surface brightness and sizes can be accommodated. 
We have one dwarf and several LSB galaxies in our sample. At the same time the model, 
based as it is on a tractable physical principle consistent with laws of mechanics and special 
relativity, does not suffer from the arbitrariness and puzzling inconsistencies of MOND. 

In order to describe a universal galaxy rotation curve [14, 15] one at most needs three 
parameters-to specify the initial slope, where it bends, and the final slope. Any more pa- 
rameter is redundant. In this regard the string model utilizes just the right number of them. 
The fact that it fits well on par with the dark matter model which employs two extra param- 
eters should be taken seriously. However one should guard against reading too much into the 
game of fitting. For example, one cannot obtain a unique decomposition of the mass compo- 
nents of a galaxy using the rotation curve data alone, a difficulty having been encountered 
in the context of comparing different dark matter halo profiles. Acceptable fits (defined as 
+ 1 [16]) can be gotten with dark matter alone without any stellar matter in CDM 
model. And the roles of dark matter and stellar matter can be completely reverted in the 
fitting routine. The physical difference is, on the other hand, dramatic. This degeneracy 
is less severe in the string model in the sense that string field cannot be completely traded 
off in favour of stellar matter, or vice versa. However there still exists a range of values for 
Rd, po and H, where "acceptable" fits can be obtained. Therefore given the quality of the 
presently available data rotation curve fitting alone cannot distinguish between dark matter 
and string field in galaxies. 

However precision measurements extended to radii r ~ 20R(i can distinguish string 
model from the other models: a gently rising rotation curve in this region is a signature 
prediction of this string toy model. At this moment we are, nevertheless, encouraged by this 
inchoate results to pursue further. In a separate article we shall subject our string model 
to other reality checks, and we shall report on how this simple string model accounts for 
gravitational lensing which is often quoted as another strong evidence for the existence of 
dark matter at intergalactic scales. 

At this point it is worth mentioning that a critical reanalysis of available data on velocity 
dispersion of F-dwarfs and K-giants in the solar neighborhood, with more plausible models, 
performed by Kuijken and Gilmore concluded that the data provided no robust evidence for 
the existence of any missing mass associated with the galactic disc in the neighborhood of the 
Sun [17]. Instead a local volume density of po = 0.10MsunPC~^ is favored, which agrees with 
the value obtained by star counting. Dark matter would have to exist outside the galactic 
disk in the form of a gigantic halo. Their pioneer work was later corroborated by [18-21] 
using other sets of A-star, F-star and G-giant data. Note that this observation can be nicely 
explained by our model as the field only affects the centripetal motion on the galactic plane 
but does not affect the motion perpendicular to the galactic plane. 

We presented a simple string toy model with only one free parameter and we showed 
that it can fit the galaxy rotation curves equally well as the dark matter model with the 
generalized NFW profile. The latter employs two more free parameters compared with the 
string model. The string model respects all known principles of physics and can be derived 
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from first principle using string theory, which in turn unifies gravity with other interactions. 
Our model has an unambiguous prediction concerning the rotation dynamics of satellites and 
stars far away from the center of the (host) galaxy. The ability to test the validity of string 
theory as a description of low energy physics makes the exercise worthwhile. 

2 Perihelion precession 

In this section we test the string model with planetary precession data in the solar system. 
In [22], anomalous precessions for Mars and Saturn was reported, for which no explanation 
within the the standard paradigm seems to exist. It is thus an interesting laboratory to 
test our string model and ask if the string gauge field can explain the reported anomalous 
precession. Furthermore, this also serves as a check of the string model independent of 
the galaxy rotation curves. It is because we can use the data of anomalous precession to 
determine the profile as well as the strength of the string gauge field inside the solar system. 
As it turns out, the extra field needed to generate the extra centripetal force to account for 
the anomalous precession has a profile of a dipole field generated by the Sun. It is then 
important to compare the magnitudes of field strength as obtained by different methods and 
observations. A consistent model should give similar values in the field strength for the same 
object. 

2.1 Field strength from precession in solar system 

Here we will use the anomalous precession data to determine the string field profile and field 
strength in the solar system. We attribute all the anomalous precession to the magnetic like 
force due to the string field. We are interested in the quantity [2] 

n = 9E (2.1) 

m 

where ^ and H are the string charge-to-mass ratio and the field strength, respectively. ^ 
Q has dimension, s~^, that of frequency. In other words we are testing the validity of 
Newtonian gravity in the extremely low frequency regime. The corresponding quantity in 
electromagnetism is — . In these units it is easy for us to compare it with strengths of other 
magnetic-like forces, e.g in electromagnetism. Precession from a magnetic-like force 
perturbation has been worked out in detail from first principle in [23] 




To calculate the strength of Q using precession data we replace ^ with Q = in the above 
formula (2.2) and invert it to get 

n = {-5Cj)^GMq^ (2.3) 
7ra2 

Relative errors from other contributions are extremely small compared to the one of the 
precession rate, thus the error of can be computed by 

Err{n) = Err{6u)^/GM^^ , (2.4) 



^Note that here we dropped the ^ factor in Q as defined when discussing GRC fitting as constants of order 
one are not important. 
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Figure 6. $1(1 x 10 s ) versus distance r{km) to the Sun 



or simply 



Err{n) 




Err{5uj) 


n 




5Co 



(2.5) 



The value of Vt determined from precession data are shown in the following table, Table 1. And 
the fitting results are shown in Fig. 6. Orbital data used in fitting are from HORIZON [24] . 
Precession data are obtained from [22]. 

The central values of for the inner planets exhibit a decreasing pattern with respect 
to r, although big error bars allow also for the case of vanishing string field. At Saturn, Q. 
is nonzero within one a. However, as mentioned in [22], the error bar at Saturn may actually 
be bigger than quoted, in which case the value of precession, or may vanish. More precise 
measurements on precession are needed to determine definitely the existence of (in other 
words, anomalous precession rate) in the solar system. Inspired by this decreasing pattern in 
the section 2.2 we will fit data of inner planets with a (nearly) power law profile. No matter 
how critically we take the profile and magnitude of here, nevertheless it is certain that the 
upper limit of O at the solar system is on the order of 10~^'^s~^. 

As a comparison, let us note that for the real magnetic field near earth, B ~ 10~^ Tesla, 
and with ^ = 1.76 x lO^^C/kg, we have flem ~ 100s~^. (See Appendix A for details.) In 
that sense, the string field strength is, naively, 10^^ times smaller than the strength of the 
magnetic field near the earth, assuming the charge to mass ratio under the string field is 
the same as that of an electron. We may wonder why this "strong" magnetic field has not 
affected precession of planets, specifically we can ask if it is related to the observed anomalous 
precession of planets. One reason why we do not have to worry about the real magnetic field is 
that planets are electrically neutral, but charged under the string gauge field as postulated. 
The real magnetic field can act on neutral matter only through dipole-dipole interaction, 
which, as explained in the Appendix A, for several reasons, does not contribute to planet 
precession. 



The same argument applies to precession rate as well, which is proportional to il. 



Name 


Mercury 


Venus 


Earth 


Mars 


Saturn 


<5cj(10-4 "/cy) 


-36 ± 50 


-4 ±5 


-2 ±4 


1 ±5 


-60 ± 20 


T(y) 


0.240846 


0.6151970 


1.0000175 


1.8808 


29.4571 


a(km) 


57909100 


108942109 


152097071 


249209300 


1513325783 


e 


0.205630 


0.0068 


0.016710219 


0.093315 


0.055723219 


Jl(lO-i's-i) 


1.11 


1.22 X 10-1 


5.99 X 10-2 


-2.69 X 10-2 


1.68794 



Table 1. : Planets data and values of f2 



There is however another question we may ask about the string model: now we are 
assuming matter is charged under string field and being acted on by the corresponding 
magnetic part for this charge, then is there an electrical part of interaction between stringly 
charged matter? After all, so far we seem to be assuming all matter take the "same" kind 
of charge. This question lies outside of our current model and calls for more theoretical 
investigations into the nature of this string charge. As for the model used here, we are 
considering a magnetic interaction in the form of a Lorentz force. ^ 

A few words on the Saturn are warranted. One aspect special about Saturn is that it 
belongs to the gas giant group while all other planets considered here are small and solid 
planets of the inner planet family of the solar system. As in the electromagnetic theory, the 
content and structures of planets may affect their interactions with the string field, which 
might explains the somewhat anomalous behavior of the Saturn. This speculation could be 
supported if anomalous precession behaviors of other gas giants similar to the Saturn can be 
measured in the future. 

2.2 Profile fitting of precession in solar system 

The decreasing pattern of for inner planets with respect to r indicates it might be useful 
to fit these values with a power law term. Presumably we could attribute this r dependent 
term to the Sun from which r is measured. On the other hand, note that at the Mars 
is negative, although it seems also sitting on the same curve passing the first three inner 
planets. 

One possible configuration for Vt compatible with these 2 facts is then that, in addition 
to a power law term, there is also a weak constant background with opposite direction to 
the from the Sun. The constant background might be provided by all other matter in the 
universe. The Milky Way should be the most important source of this influence. Thus here 
we try to fit field strength at different inner planets with the following profile, 

9.{r) = A + Br-" (2.6) 

For the actual fitting on computer, the profile used is 

^="» + "'(lok)"° 

^ It is amusing to find out the extra acceleration due to the string field for objects on Earth. Since the field 
strength is gotten for the rest frame relative to the Sun, the velocity of objects on the earth should be nearly 
the same with the velocity of Earth relative to Sun, i.e 30fcm/s. The corresponding acceleration produced is 
therefore ~ lQ~^^m ■ . 
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Planet 


Mercury 


Venus 


Earth 


Mars 


ratio 


50.3 


7.1 


2.5 


0.54 



Table 2. relative strength: power law term/constant term 



The best fitting parameters for this profile is 

no = -0.0223787 x lO^^^^s^^ (2.8) 
= 260.504 X lO^^'^s"^ 
a = 3.09953 

It is helpful to know the relative strength of the two components from the Sun and the 
constant background, which is shown in Table 2. The power law term decreases with r 
quickly relative to the background. We can safely say that in most areas in the solar system, 
the string field would just be around that background, which is on the order of 10~^^Hz. 
And that a is found to be near 3, which is exactly the power for a dipole field. It means that 
the string field interaction between the sun and planets is similar to that between a magnetic 
dipole and charged particles. 

Also note that the fitting result tells us that the constant background in the solar system 
is negative. This is good news for the string model. As in electromagnetism, we expect the 
string gauge field in a galaxy to be generated by the rotation of (stringly charged) matter in 
the galaxy, just like rotating electric charge would generate a magnetic field. Since the Sun 
and planets in the solar system rotate in the opposite direction of that of stars' rotation in 
the Milky Way, it is therefore reasonable to expect the background field to be negative if we 
consider the one from the Sun as positive. This is exactly what the profile fitting has told 
us! Here only data of the solar system was used, but the conclusion is for the whole Milky 
Way, specifically for its rotation direction. 

However, this conclusion should be taken with a grain of salt, as we will explain below. 
Firstly, the long error bars of precession weaken the conclusion from this profile fitting. 
Secondly, there exists possibilities for the field from the Sun to be in the same direction with 
that of the Milky way. As in ordinary electromagnetic theory, for a right handed current disk 
the magnetic field is downward outside of the major current distribution, but upward if we go 
into the current disk somewhere, specifically at the center of the disk. There is a place within 
the concentration where the magnetic field changes its sign. ^ Thus even if the Milky Way is 
rotating in opposite direction from the solar system, if the Sun is too close inside the major 
mass concentration of the Milky Way, the background therefore should still be positive. To 
our advantage, it is known that the Sun contains 99% of the total mass in the solar system, 
so it is reasonable to assume all planets are well outside of most mass in the solar system. 
And the solar system lies somewhat outside the majority of mass concentration of the Milky 
Way. 

2.3 Field strength from Milky Way rotation curve 

In section 2.1 we obtained Q in the solar system from anomalous precession. The constant 
background value, found by the profile fitting in section 2.2, should mainly comes from other 
matter in the Milky Way altogether. On the other hand, by the same idea used in [25], we 

^Considering this, the f2 in [2, 25] are position-averaged one over the galaxy. 
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can also estimate Q at the Solar system due into the Milky way's rotation. Thus a natural 
check of the string model is to compare these two values of O: one of the constant background 
from precession in the solar system, another one from the rotation curve of the Milky Way. 

Directly using the idea in [25], we can make a rough estimate of in the Milky way 
as follows. In the string model, the total force on the galaxy mass is composed of only the 
gravitational attraction from visible mass in the Milky Way and the magnetic-like force from 
the string field. With increasing r, the gravitational force decreases quickly, and the magnetic 
like force always increases. At the position of the Sun the rotation curve is well into the flat 
region [26-28] . Therefore on the Sun the gravitational force should be negligible with respect 
to the magnetic-like force from the string field. Then [2] 

m— ^ QHv (2.9) 
r 

and thus 

n^- (2.10) 

r 

For the sun [26-28], Vq ~ 200/cm • and tq ~ 7.Qkpc. So 

w 26s"^ • ^ = 8.5 X lO^^^s"^ (2.11) 
kpc 

This is the upper limit on $7 due to the Milky way at the position of the Sun. Since 
there is still a portion of distance further out to nearly 20kpc where the curve is quite flat 
(with velocity ~ 200km • s~^)[26-28], we could have used these distances instead of 7.6kpc 
and the value of $7 will be reduced by a factor of about three. In any case, it is safe to say 
the upper limit of Q is on the order 10~^^s~^. This is the strength of the field component 
perpendicular to the galactic plane. To convert it to the solar system, notice that the north 
galactic pole and the north ecliptic pole form an angle of 60.2° (which means the field in 
the solar system would be reduced almost by half), and the milky way rotates clockwise 
when viewed from north galactic pole (which means the field is negative in the solar system 
since planets in the solar system rotate counterclockwise if viewed from the same direction). 
Therefore from the rotation curve of the milky way the upper limit of the effective field 
strength in the solar system due to matter in the milky way is — 10~^'' ~ — 10~^®s~^. 
This magnitude is close to the one found by direct precession calculation without the profile 
assumption, but it is not the case for field direction. The precession indicates the field in 
the solar system is positive, while rotation curve of the milky way says it is negative. One 
possible explanation is that at places near the Sun the field is dominated by the positive field 
generated by the Sun, and the Milky way provides only the constant background in the solar 
system. The profile fitting of precession, apart from a dipole like part due to the Sun, indeed 
gives a negative background (—0.02 x 10~^''s~^). However the magnitude there is smaller by 
almost two orders of magnitude. 

2.4 Discussion 

With anomalous precession data, we calculated the upper limit on the field strength at several 
places in the solar system. It is found to be on the order of 10^^'^s^^. The distribution of 
the field in the solar system looks like a superposition of a constant background and a 
decreasing "dipole" component in (2.6). We discussed a possible configuration that the Milky 
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way produced the background and the Sun sources the "dipole" component. A profile fitting 
is done for this configuration. The background is found to be negative, which, when combined 
with analysis using electromagnetism analogy, correctly matches with the fact that the solar 
system and the Milky Way rotate in opposite directions (assuming the solar system is far 
enough away enough from the Milky way center). 

As a comparison with the field strength gotten from precession in the solar system, we 
estimated the upper limit on the field strength of the string gauge field in the Milky way 
directly from rotation curve of the Milky Way at the position of the Sun. This estimate is 
of a similar order of magnitude with the one gotten from precession. But the direction is 
opposite. If we consider the background plus dipole field configuration used for profile fitting 
to be correct, then we should consider only the background value when comparing with the 
one estimated by Milky way rotation curve. In that case, although the background direction 
agrees with the one predicted by Milky Way rotation curve, but the magnitude is smaller by 

2 orders of magnitude. 

3 Summary 

We discussed the possibility that the rotational speed of the stars about the center of a spiral 
galaxy is supported by the presence of string gauge field which couples universally to all forms 
of matter. We compared the good of fit of the string model with the commonly accepted 
Cold Dark Matter model with the generalized Navarro-Frenk- White profile for Dark Matter 
distribution. We fit rotational speed data of 22 spiral galaxies of varied range in size 
and luminosity. Dark Matter model fit marginally, 9%, better at the price of two more free 
parameter than the string model. A Tully-Fisher relation relating the visible mass and field 
strength of the string gauge field to velocity GMstar^ = can be derived dynamically from 
the string model, which is obeyed fairly well by all 22 galaxies of varied sizes and luminosities. 

The existence of string gauge field is taken a step further and is applied to explain the 
anomalous precession of planets in the solar system. The extra field needed to generate the 
extra centripetal force to account for the anomalous precession has a profile of a dipole field 
generated by the Sun. 

The values of string fields in the following cases 

• galaxy rotation curves of a set of 22 galaxies (without the Milky way), 

• perihelion precession of some planets in the solar system and 

• Milky way rotation curve 

are summarized in Table 3. Interestingly, and also luckily for the string model, results from 
these different methods lead to similar order of magnitude for 0,. 

Acknowledgments 

We would like to thank Yuran Chen and Youhua Xu for collaboration at the early stage of the 
precession project. We also thank Konstantine Savvidy for helpful discussions related to the 
rotation curve project. We have also benefited from discussions with Gang Chen, Zheyong 

^°The NFW profile with a = 1 is Icnown to have difficulty fitting dwarf galaxies as well as galaxies with low 
surface brightness. Hence the generalized profile which leaves a a free parameter is called for. 



-17- 





Solar System C Milky Way 


22 other galaxies 


Rotation Curve 


\n\ < io-i« 


6 X IQ-^^ < \n\ < 10"^^ 


Precession 







Table 3. Dimension for O is s ^. The first column indicates the kind of observation used, and the 
first row indicates the object on the which the observation is done. 
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A Solar system magnetic field and its eff"ect on planet perihelion precession 

References for this appendix are [29-33]. The Sun and most planets in the solar system have 
magnetic field due to dynamo effect. If we treat both the Sun and the planet as magnetic 
dipoles interacting in vacuum (which leads to a central force with n = —4), using data of 
magnetic fields of the Sun (around 1 ~ 2 gauss at the polar region) and the Earth (around 
0.6 gauss at the polar region), we can find the corresponding precession produced is nearly 
10~^arcsec/cy, which is 2 orders of magnitude smaller than the observed one. In fact the 
magnetic field in the solar system is much more complicated than those produced by several 
dipoles in vacuum. First, the solar system is not empty but filled with particles emitted from 
the Sun, i.e the solar wind. Charged particles lock with it the magnetic field of the Sun and 
spread it all around in the solar system. From the Sun to about the position of the Earth, 
the magnetic force line is parallel to the radial stream of solar wind particle and falls off 
by r~^. From the position of the Earth to about position of the Mars is a field free region 
with B < 10~^ gauss. Further out to the position of the Jupiter is a region with disordered 
magnetic field with B ~ 10~^ gauss. For precession, the most important feature of the solar 
system field is that it is oscillating. First, for the Sun the magnetic dipole axis is inclined 
relative to the rotational axis. This leads to an oscillating neutral current sheet. Therefore 
planets on the ecliptic plane is above the neutral current sheet for half of solar self rotation 
period, below for another half. Since field directions above and below the neutral current 
sheet is opposite, the magnetic force experienced by the planets also change directions within 
one self rotation of the Sun. Second, the magnetic field of the Sun also changes direction every 
22 years due to its differential rotation, which leads to another oscillation of magnetic field 
on the planets. Altogether these two oscillations make the magnetic field effect on precession 
neglectable with respect to other accumulating effects. 

B Different fitting models give difi'erent best fittings 

Fitting a galaxy consists of following steps: 

• assume a density profile (a parametric description of the density), 

• adjust the values of free parameters to minimize an "error function." 

The corresponding resulting parameters are called best fitting parameters. For a single 
galaxy we can define different error functions. It can be defined by rotation curve, by surface 
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brightness or some other observation data. The point of this appendix is that we usuahy 
get different best fitting parameters when using different error functions. In particular, best 
fitting parameters for rotation curve are different from those for surface brightness. Below 
we provide a simple and idealized example to illustrate this point. 

Real Suppose the real density distribution is a linearly decreasing function of radius and 
becomes zero outside of a cutoff radius, 

where po and R are two fixed constants for this particular galaxy. The corresponding 
velocity, using Newtonian gravitation theory, is 



j ^2npoGM{'^-lj^)r,iO<r<R), 

v{r) = < ^ / 

[ ^27TpoGNlf, ir>R) 

and brightness is (assuming light is proportional to mass) where 7 is light to mass ratio. 
Which gravity theory we use does not affect the conclusion of this appendix, so long as we 
use the same theory for any profile to derive the velocity. 

Guessed Without knowing the real distribution, suppose we assumed for this galaxy a 
constant distribution profile 

Here a and D are two parameters, rather than constants, being fitted to get best fitting 
values, while po and R have particular fixed values for this galaxy. This density profile 
produces following velocity profile 

r V^aG^r, (0 < r < D), 
"'"-^'[f^^,, (r>D). P'-^) 

and brightness profile 

- i 0, (r > D). (^-^^ 

Fitting As said before, we can define different error functions to do the fitting. We can use 
either rotation velocity or brightness for fitting. Ideally the error functions can be defined 
respectively for rotation velocity and surface brightness by 

POO 

xl{D,a)= dr'[vmodeiir',D,a) - Vreaiir')f (B.7) 
Jo 

POO 

xl{D,a)= / dr'[Brnodei{r,D,a)-Breaiir')f (B.8) 
Jo 



^'^Here we assume the galaxy is a disk and the distribution has only r dependence. So it is actually more 
appropriate to call it surface density. 
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After minimizing these two error functions we get two sets of best fitting parameters for D 
and a. Let's denote the best fitting value for velocity error function by {Da,cfa)-, for 
brightness error function by {Db,(Tb)- Below we show these two sets of values do not 
coincide. 

If we use the error function for brightness, obviously the best fitting parameters are 

Db = R, <Jb = ^Po (B.9) 

To show (DbjCtb) 7^ {Da,(^a) it is sufficient to show that {DbjCtb) does not minimize xt- 
For simplicity we take units such that 

po = l, R = l, 7 = 1, 7rGjv = l, (B.IO) 

then 

Db = 1, (TB = l- (B.ll) 



2 



In this unit system the real velocity is 



%, (0<r<l), 

(^>i). 

Taking 100000 as the upper limit of the integration, we find 

X^.pB = l,cTB = 0.5) = 0.2, (B.13) 
xl{D = 0.9, a = 0.5) = 0.05 (B.14) 

Thus {DbiCtb) 7^ {Da,(ta)- The best fitting parameters for brightness do not best fit the 
velocity curve. 

Remarks Above wc considered a simple and idealized galaxy and showed that best fitting 
parameters for different error functions are in general different, although we were doing the 
fitting for the same galaxy. Another point worth noting here is that if we had guessed at the 
correct density profile (which linearly decreases and vanishes beyond some cutoff radius) , the 
two best fitting parameters will be the same. Wc get different fitting results because we used 
a "wrong" profile for this galaxy. In real life the mass distribution for a galaxy is extremely 
complicated and can not be exactly described by any simple "profile function." Hence after 
we assume the profile, define the error function and then do the fitting, we will always get 
different fitting results for "independent" error functions (e.g velocity and brightness). If the 
guessed profile is closer to the real distribution we get closer results for fittings by different 
error functions. In other words, a big difference between fitting results by different error 
functions means the profile we guessed at is very different from the real one. 

Given different profile assumptions for the galaxy distribution, we thus have a way to 
judge in some sense which one is more "correct". With each profile we can derive corre- 
sponding distributions of velocity, brightness and so on, and with each of these distributions 
we can compute( if we have the data) an error function x?, which is dependent on a set of 
parameters owned by this profile, 

, Velocity xl 
Density Profile — > { Brightness x| 
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galaxy measurement 







name 
mSldwb 
ngc0628 
ngc0925 
ngc2403 
ngc2976 
ngc3031 
ngc3184 
ngo3198 
ngc3521 
ngc3621 
ngc3938 
ngc4236 
ngc4321 
ngc4536 
ngo4569 
ngo4579 
ngc4625 
ngc4725 
ngc5055 
ngc5194 
ngc6946 
ngc7331 



B magnitude Distance(Mpc) 



mucin 

N/A 
29.95 
29.85 
27.68 

28.1 
28.57 

30.2 
30.42 
30.29 
29.58 
30.81 
27.08 

31.9 
32.1 1 
30.52 

31.8 
30.35 
31.45 
30.09 
30.01 
29.12 
30.75 



2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 



-12.5 3.5 

-20.60 9.77 

-20.05 9.33 

-19.56 3.44 

-18.12 4.17 

-21.54 5.18 

-19.88 10.96 

-20.44 12.13 

-21.08 11.43 

-20.51 8.24 

-20.01 14.52 

-18.10 2.61 

-22.06 23.99 

-21.79 26.42 

-21.10 12.71 

-21.68 22.91 

-17.63 11.75 

-21.76 19.50 

-21.20 10.42 

-20.51 10.05 

-20.89 6.67 

-21.58 14.13 



Figure 7. Relevant galaxy observation data for determining luminosity and distance. 

If the profile perfectly match the real one, there exists a single set of parameters simultane- 
ously making all \^ zero. On the other hand, if the profile differs too much from the real one, 
even if we can make one of the small, the corresponding parameters will usually make 
other x^'s very large. Therefore it makes sense to use, e.g. the average of several x^'s as 
the error function to be minimized. This profile, with the corresponding minimizing param- 
eters of the average error function, best describes the distribution averagely, i.e. considering 
distributions whose x^ is averaged. 

However we cannot use this to find the real distribution. We can only compare profiles 
given their profile assumptions and say which is better in describing velocity, brightness and 
so on, or if we use some averaged error function, also say which is better considering their 
general performances in describing various properties simultaneously. 

C Detailed Fitting Results of 22 galaxies using dark matter model and 
using string model 

In this appendix we present all the results on the data fitting. Figure 7 summarizes the 
relevant data for luminosity as well as for distance determination. In Figure C and Figure 9 
the best fit values for the five free parameters in the dark matter fitting and three free 
parameter in the string model fitting are presented. Mass of the stellar mass and dark matter 
halos (in the case of dark matter model) as well as the mass to light ratios are computed from 
the best fit values for each galaxy. These values are tabulated in Figure C (dark matter) and 
in Figure 11 (string). The rest of the appendix presents 22 graphs of dark matter fit and 
string fit side by side for each of the 22 galaxies. In each of the graph small cubes with error 
bars represent observational data the other symbols, and the curves, represent theoretical 
values. The X-axis is radius in k-pc and Y-axis is velocity in kms~^ . 
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Figure 8. Best fit values for the five free parameters in the Dark Matter model with generalized 
NFW profile. 
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Figure 9. Best fit values for the tliree free parameters in tire 


string model. 
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Figure 10. Summary of stellar mass and mass of dark matter halos computed from the best fit 
values, as well as mass to light ratios of the 22 galaxies in the dark matter model. The mass is in 
units of Mgun- Mass to light ratio is relative to the mass of light ratio of the Sun. 
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Figure 11. Summary of stellar mass as well as mass to light ratios as computed from the best fit 
values in string model. The mass of the galaxies is in units of Ms„„. Mass to light ratio is relative to 
the mass of light ratio of the Sun. 
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Figure 12. Rotation curve of dwarf galaxy mSldwb fit with the dark matter model (left) and with 
the string model (right). The X"Squared value per degree of freedom using the dark matter model 
with a NFW profile is 0.158 while that using the string model is 0.1254. 
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Figure 13. Rotation curve of ngc0628 fit with the dark matter model (left) and with the string model 
(right). The ^-squared value per degree of freedom using the dark matter model with a NFW profile 
is 4.398 while that using the string model is 4.3468. 
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Figure 14. Rotation curve of ngc0925 fit with the dark matter model (left) and with the string model 
(right). The x-squared value per degree of freedom using the dark matter model with a NFW profile 
is 2.392 while that using the string model is 2.4518. 




Figure 15. Rotation curve of NGC2403 fit with the dark matter model (left) and with the string 
model (right). The x-squared value per degree of freedom using the dark matter model with a NFW 
profile is 4.467 while that using the string model is 4.2468. 
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Figure 16. Rotation curve of ngc2976 fit with, the dark matter model (left) and with the string 
model (right). The ^-squared value per degree of freedom using the dark matter model with a NFW 
profile is 0.504 while that using the string model is 0.4001. 
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Figure 17. Rotation curve of ngc3031 fit with the dark matter model (left) and with the string 
model (right). The ^-squared value per degree of freedom using the dark matter model with a NFW 
profile is 5.614 while that using the string model is 5.6977. 




Figure 18. Rotation curve of NGC3184 fit with the dark matter model (left) and with the string 
model (right). The ^-squared value per degree of freedom using the dark matter model with a NFW 
profile is 0.573 while that using the string model is 0.5740. 
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Figure 19. Rotation curve of ngc3198 fit with the dark matter model (left) and with the string model 
(right). The x-squared value per degree of freedom using the dark matter model with a NFW profile 
is 0.441 while that using the string model is 0.2748. 




Figure 20. Rotation curve of ngc3521 fit with the dark matter model (left) and with the string model 
(right). The x-squarcd value per degree of freedom using the dark matter model with a NFW profile 
is 0.167 while that using the string model is 0.5875. 




Figure 21. Rotation curve of NGC3621 fit with the dark matter model (left) and with the string 
model (right). The x-squarcd value per degree of freedom using the dark matter model with a NFW 
profile is 0.125 while that using the string model is 0.3658. 
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Figure 22. Rotation curve of ngc3938 fit with tiie darlc matter model (left) and witli the string model 
(right). The x-squared value per degree of freedom using the dark matter model with a NFW profile 
is 0.905 while that using the string model is 1.0304. 
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Figure 23. Rotation curve of ngc4236 fit with the dark matter model (left) and with the string model 
(right). The x-squarcd value per degree of freedom using the dark matter model with a NFW profile 
is 1.180 while that using the string model is 0.3255. 




Figure 24. Rotation curve of NGC4321 fit with the dark matter model (left) and with the string 
model (right). The x-squarcd value per degree of freedom using the dark matter model with a NFW 
profile is 0.956 while that using the string model is 2.0703. 
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Figure 25. Rotation curve of NGC4536 fit witli tlie dark matter model (left) and with the string 
model (right). The ^-squared value per degree of freedom using the dark matter model with a NFW 
profile is 0.598 while that using the string model is 0.7394. 




Figure 26. Rotation curve of ngc4569 fit with the dark matter model (left) and with the string model 
(right). The x-squared value per degree of freedom using the dark matter model with a NFW profile 
is 0.651 while that using the string model is 0.6621. 





Figure 27. Rotation curve of NGC4579 fit with the dark matter model (left) and with the string 
model (right). The ^-squared value per degree of freedom using the dark matter model with a NFW 
profile is 0.719 while that using the string model is 0.688. 
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Figure 28. Rotation curve of ngc4625 fit with, tlie dark matter model (left) and with the string 
model (right). The ^-squared value per degree of freedom using the dark matter model with a NFW 
profile is 0.874 while that using the string model is 0.5478. 




Figure 29. Rotation curve of ngc4725 fit with the dark matter model (left) and with the string model 
(right). The x-squared value per degree of freedom using the dark matter model with a NFW profile 
is 0.481 while that using the string model is 0.2683. 




Figure 30. Rotation curve of NGC5055 fit with the dark matter model (left) and with the string 
model (right). The x-squared value per degree of freedom using the dark matter model with a NFW 
profile is 1.133 while that using the string model is 3.5208. 





Figure 31. Rotation curve of ngc5194 fit with the dark matter model (left) and with the string model 
(right). The x-squared value per degree of freedom using the dark matter model with a NFW profile 
is 0.442 while that using the string model is 0.8067. 




Figure 32. Rotation curve of ngc6946 fit with the dark matter model (left) and with the string model 
(right). The x-squared value per degree of freedom using the dark matter model with a NFW profile 
is 5.670 while that using the string model is 6.0346. 
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Figure 33. Rotation curve of NGC7331 fit with the dark matter model (left) and with the string 
model (right). The x-squared value per degree of freedom using the dark matter model with a NFW 
profile is 0.839 while that using the string model is 0.5221. 
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